Casimir scaling in G2 lattice gauge theory 
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We computed potentials between static color sources from the six lowest representations of G2 
lattice gauge theory, in numerical simulations with the Wilson action on asymmetric lattices with 
nonperturbatively estimated values of the bare anisotropy. We present evidence for (approximate) 
Casimir scaling of the obtained intermediate string tensions. The agreement with the Casimir-scaling 
prediction improves by increasing the coupling j3 in the weak-coupling region above the crossover 
observed in G2 gauge theory. The result naturally fits into confinement models with magnetic 
disorder and vacuum domain structure, but may represent a challenge for other approaches. 
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I. INTRODUCTION 

Any theory aiming to explain the mechanism of color 
confinement in quantum chromodynamics has to face 
a would-be simple - but rather non-trivial in reality - 
task: to describe qualitative features of potentials be- 
tween static charges from various representations of the 
gauge group. If one neglects dynamical quarks in the 
first approximation, the static quark-antiquark poten- 
tial in SU(iV) gauge theory exhibits distinct behavior 
in three ranges of intcrquark distances. At short dis- 
tances, the interaction is dominated by gluon exchange 
and the potential is Coulomb-like, its strength being pro- 
portional to the eigenvalue of the quadratic Casimir op- 
erator in the given representation of color sources. At 
intermediate distance scales, from the onset of confine- 
ment to the onset of color screening, the potential is ex- 
pected to be linearly rising, and the corresponding string 
tensions for different representations are again approxi- 
mately proportional to the quadratic Casimir. 1 This ef- 
fect, dubbed Casimir scaling in Ref. 0L was observed in 
numerical simulations of both SU(2) [!, H, Q and SU(3) 
@; H, B EH EH lattice gauge theories. 2 Finally, at asymp- 
totic distances, the string tensions depend only on the N- 
ality k of the representation: representations with k = 
are screened and their potentials are asymptotically flat, 
k =/= representations are confined and their potentials 
continue to rise linearly. 

A confinement scenario which naturally accounts for 
both Casimir scaling at intermediate distances and the 
asymptotic TV-ality dependence is based on condensation 
of center vortices in the QCD vacuum (see [l2| for a re- 
view and summary of numerical evidence for the model). 
The vacuum is filled with percolating thick vortices car- 
rying color magnetic flux quantized in terms of elements 
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1 This proportionality is exact in two dimensions, and can be ex- 
tended to D = 3 and 4 by dimensional reduction arguments, see 
e.g. QUEi- 

2 Ref. Ildl contains an extensive bibliography of earlier references. 



of the gauge group center. Asymptotically large Wilson 
loops, linked to center vortices, "feel" the total (center) 
flux; this is how iV-ality shows up easily in asymptotic 
string tensions. On the other hand, Wilson loops of in- 
termediate size cross thick-vortex cores and enclose only 
a part of their magnetic flux. As a consequence, inter- 
mediate string tensions can exhibit approximate Casimir 
scaling. A simple model of how and why it happens, 
was suggested in Ref. (lH ] (some numerical support was 
provided by flij). 

Recently, G2 Yang-Mills theory has attracted consid- 
erable attention [TEEl, EE El EE SI El, H as a means 
to look at the physics of confinement from a fresh new 
angle. The G2 group has a trivial center: the center of 
its universal covering group is also trivial [15j, and there 
seems to be no reason to believe that the group center 
plays any important role in G2 gauge theory at all. In 
the strict sense of the word the theory is only temporar- 
ily confining, since potentials for static color charges from 
any representation, including the fundamental one, must 
be asymptotically flat. Unlike SU(AT), even a "quark" 
from the fundamental ({7}) representation of G2 can be 
screened by a bunch of "gluons" (from the 14-dimcnsional 
adjoint representation). However, this fact does not re- 
ally contradict the center-vortex confinement scenario: 
G2 does not possess non-trivial center vortices, and the 
asymptotic string tension is zero if they are not present. 
As argued by the authors of Ref. [2l[ , the G2 temporary 
vs. permanent confinement issue is more semantics than 
physics; the really important point is that even in the 
G2 gauge theory one expects the static potentials grow 
linearly over a certain range of distances, from the scale 
where perturbation theory breaks to the onset of screen- 
ing. The linear rise of the fundamental potential in this 
intermediate region was clearly demonstrated in numer- 
ical simulations [2l[ ■ 

All asymptotic string tensions are zero in G2 - what 
about Casimir scaling at intermediate distances? Dimen- 
sional reduction could be invoked to argue in its favor. 
However, if it holds true and if finite thickness of center 
vortices was able to explain successfully Casimir scaling 
in SU(iV), what can it be attributed to in a theory with 
trivial center? Is there any origin of Casimir scaling com- 
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mon for "centerless" and "centerful" gauge models? The 
answer given in [2~fl | was affirmative, the authors pro- 
vided a model - an extension and improvement of that 
of Ref. [HI - in which the Yang-Mills vacuum state has 
a domain structure, with the color magnetic flux in each 
domain quantized in units of the gauge group center (be 
it trivial or not), and Casimir scaling results from random 
spatial fluctuations of the flux in each domain. 

The aim of the present paper is to subject the Casimir- 
scaling hypothesis and predictions of the model of [2l[ to 
a test in numerical simulations of G2 gauge theory on a 
four-dimensional Euclidean lattice. First, in Sec. [Til wc 
will briefly summarize elements of the model. Then, in 
Sec. IIII) we sum up a few facts about the exceptional G2 
group, explain the applied lattice methods, and present 
a subset of technical details on the enhancement of the 
ground-state overlap via smearing, and the determina- 
tion of renormalizcd anisotropics. Section HVl summarizcs 
results, in particular subsection IIV Bl the fundamental 
representation potential, while HVCl potentials for higher 
representations. We conclude with a brief discussion and 
summary (Sec. [V}. A few technicalities are relegated to 
appendices. 



II. A MODEL OF THE YANG-MILLS VACUUM: 
RECAPITULATION 

An essential element of the derivation of Casimir scal- 
ing in the D = 2 Yang-Mills theory is the following: the 
YM action is quadratic in the field strength tensor i 7 ^, 
contains no derivatives of or constraints on Fij, there- 
fore field strengths of vacuum fluctuations in different 
points are uncorrelated, and disorder all-representation 
Wilson loops. This cannot be true for an effective action 
in D > 2 dimensions. Even though the leading term of 
a derivative expansion of the Yang-Mills effective action 
is again quadratic in field strengths [l|, 3 the fields are 
constrained by Bianchi identities and higher-order cor- 
rections contain derivatives. This is not a drawback, but 
desired: constraints and derivative terms must somehow 
contrive, in D > 2 dimensions, to change the propor- 
tionality of string tensions to quadratic Casimirs to the 
asymptotic Y-ality dependence. 

In the simple model of Ref. [2l[ it is assumed that if we 
take a 2D slice of the four-dimensional volume, we can 
split it into (partly) overlapping domains ( "patches" ) of 
a typical area Ad- Within each domain color magnetic 
fields fluctuate randomly and (almost) independently, 
with a short length of correlation I. Each domain is a 
bunch of small independently fluctuating subregions of 
area I 2 <C Ad- The only weak constraint on fluctuating 
fields bounds the total integrated magnetic flux over each 
domain to correspond to an element of the center of the 



gauge group. In SU(2), there are two types of domains, 
of the center-vortex type and of the vacuum type. The 
former correspond to the nontrivial element of the cen- 
ter, —I, and represent a cross section of a thick center 
vortex; the latter carry a zero total magnetic flux. In G2 
all domains will be of the vacuum type, since the center 
contains the identity element only. How to distinguish 
overlapping vacuum domains in 2D snapshots of the G2 
vacuum remains an unresolved task for future reflection 
and/or numerical experiments. 

Without going into details, let's briefly summarize con- 
clusions of the model. If the center of the gauge group 
contains N elements, there are N types of vacuum do- 
mains enumerated by the value k <G {0, 1, . . . , N — 1}, 
each assumed to appear with a probability fk centered 
at any given plaquette in the plane of the loop. The ef- 
fect of a domain (a 2D cross section of the kth vortex) on 
a planar Wilson loop is to multiply the loop by a group 
element 



g(a {k) ,S)=S exp 



St, 



(2.1) 



where {Hi} are generators of the Cartan subalgebra, S is 
a random element of the group, and angles af^ depend 
on the location of the vortex/domain with respect to the 
loop. If the domain is all contained within the loop area, 
then 



exp 



H 



(2.2) 



where is the kth element of the center, I is the unit 
element. If the domain is outside the loop, it has no effect 
on the Wilson loop, i.e. 



exp 



?(*) 



H = 



(2.3) 



For a Wilson loop from the representation r, the averaged 
contribution of a domain is simply 

Gr(a^) l dr = J-Xr (exp [**<*> ■ 7i\) I dr , (2.4) 

where d r is the dimension of the representation r and 1^ 
is the d r x d r unit matrix. 

One further assumes that probabilities to find domains 
of any type centered at two different plaquettes are inde- 
pendent, and that for loops smaller than the typical size 
of the domain, the r.m.s. of phases a is proportional to 
the area of the vortex contained in the interior of the loop. 
Then it is a simple exercise to show that, both in SU(AT) 
and G2, the static potential V r (R) of the representation 
r will be linearly rising for distances I < J? < V Ad, 
with a string tension approximately proportional to its 
quadratic Casimir: 



C r 



(intermediate distances). (2-5) 



See recent arguments in 12311 and references therein. 



For very large Wilson loops most vortices will be con- 
tained within the loop, the average phases are propor- 
tional to the corresponding total magnetic flux through 
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TABLE I: Group theoretical factors for G2. {D} is the rep- 
resentation, D its dimension, [Ai,A2] are its Dynkin coeffi- 
cients, and d{zj} = C{£>} /CV denotes the ratio of its quadratic 
Casimir to that of the fundamental representation. 



{D} 


{7} 


{14} 


{27} 


{64} 


{77} 


{77'} 


[Ai,A 2 ] 


[1,0] 


[0,1] 


[2,0] 


[1,1] 


[3,0] 


[0,2] 


d {D} 


1 


2 


2.3 


3.5 


4 


5 



the domain, and the prediction is, for R 3> \J Ad 

<J r ~ F[k r ,{fk}} (asymptotic distances) (2-6) 

for SU(iV), where k r is the iV-ality of the representa- 
tion r, and 

ay = (asymptotic distances) (2-7) 

for the G2 gauge group. The function F above depends 
also on the set of probabilities {/fc}; these can be tuned 
to get a desired form of fc-scaling of asymptotic string 
tensions (see [IH for a discussion of this point). 

Thus, the domain model of the YM vacuum with 
magnetic disorder [2l[, containing only a handful of ad- 
justable parameters, predicts Casimir scaling of higher 
representation potentials at intermediate distances both 
for SU(iV) and (centerless) G2 gauge theory. To verify 
the prediction for G2 , numerical simulations of the theory 
in lattice formulation are the appropriate tool. 



III. G 2 ON A LATTICE: PRELIMINARIES 

A. G 2 group 

G2 is the smallest of the 5 exceptional simple Lie 
groups, at the same time the simplest for which the uni- 
versal covering group is the group itself, and both have 
a trivial center [15j . The basic information on the group 
can be found e.g. in [TBI [i~7T]. a part is summarized here 
for the sake of completeness. The G2 group has rank 
2 and 14 generators. Its fundamental representation is 
7-dimensional. G2 is real and a subgroup of SO (7) of 
rank 3 with 21 generators, a group of 7 x 7 orthogonal 
matrices B with unit determinant. G2 matrices satisfy 7 
additional cubic constraints: 

Tabc = T def B da B eb B fC7 (3.1) 

Tl27 = 2l54 = Tig3 = T235 = T 264 = T374 = T576 = 1, 

where T is an antisymmetric tensor non-zero elements of 
which are listed above (up to permutations of indices). 

For a study of higher representations we will need ten- 
sor decompositions of various tensor products. Group 



theory provides the following relations: 

{7}® {7} = {1}® {7} ©{14} ©{27}, 

{7} ©{14} = {7} © {27} © {64}, (3.2) 

{14} ©{14} = {1}© {14} ©{27} ©{77}© {77'}, 

{7} ©{27} = {7} ©{14} ©{27} ©{64} ©{77}. 

Assume now H^ D ^(g) is a matrix corresponding to an 
element g G G2 in its D-dimcnsional irreducible repre- 
sentation. Then the relations (|3.2[) imply that traces of 
higher-representation matrices can be expressed through 
traces of the fundamental- and adjoint-representation 
matrices: 5 ' 6 

Tr © {27} = -1 - Tr B A - Tr B F + (Tr © F ) 2 , 

Tr P {64} = 1 + Tr D' 4 + Tr B A • Tr 3 F — (Tr B F ) 2 , 

Tr D {77} = -Tr D A - Tr D F - 2 Tr D A ■ Tr B F 

- (Tr D F ) 2 + (TrB F ) 3 , 
Tr D {77 ' } = Tr D A + (Tr D A ) 2 + 2 Tr D F 

+ 2 Tr B A ■ Tr B F - (Tr D F ) 3 . (3.3) 

The adjoint-representation matrix is constructed from 
the fundamental one: 

B A (g) = 2 Tr [B F (.g)t t a B F (g) t b ] . (3.4) 

The eigenvalue of the quadratic Casimir operator for 
the D-dimensional irreducible representation {D} of the 
G2 group can be computed from its Dynkin coefficients 
[Ai , A2I (see e.g. [2(|; we adopted the convention used 
in |25|). The dimension D of a representation of G2 with 
Dynkin coefficients [Ai, A2] is given by [27| : 

D = ^{l\-ilM-tlWl-il), (3.5) 

where 4 = ±(1 + Ai), £ 2 = |(4 + Ai + 3A 2 ), and £ 3 = 
5(5 + 2Ai + 3A2), while the ratio of its quadratic Casimir 
operator to that of the fundamental representation is [27l ]: 

Explicit values for the six lowest representations of G 2 
are listed in Table [U 



4 The decompositions were obtained using LiE, a clever computer 
algebra package for Lie group computations. For details, see l2- r | . 

5 We suppressed the argument g and used F and A instead of {7} 
and {14} to denote the fundamental and adjoint representation 
respectively. The symbol Tr is understood in the usual matrix 
sense, i.e. Tr DM = ^f =1 B™ ■ 

6 A careful reader might notice that, in the case of SU(3), a 
well-known decomposition {3} ® {3} = {6} © {3} leads to 
Tr D< 6 > = (Tr O f ) 2 - (Tr D F )*, while Refs. [3 used 
Tr D( 6 > = i [(Tr B F ) 2 + Tr (B F ) 2 ] . Both formulas (and similar 
pairs for other representations of SU(3)) can be shown equivalent 
using the Cayley-Hamilton theorem. 
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B. G2 lattice gauge theory 

We closely follow the lattice formulation of G2 gauge 
theory outlined in On links of a space-time lat- 

tice, one puts matrices Un(x) from the fundamental 
({7}) representation of G2. Those are represented by 
7x7 complex matrices in a parametrization suggested 
by [HI, which makes explicit the separation into the 
G2/SU(3) coset group and the SU(3) subgroup. 7 (De- 
tails of the parametrization can be found in Pepe's re- 
view [HI]-) Generators of the G2 algebra corresponding 
to this parametrization are listed in Appendix lAl 8 

We use the Wilson action: 
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J2 Co Rc Tr [P i0 (x)] 



x,i>0 



3 y, r RcTr 



(3.7) 



x,i>j>0 



where P VT {x) = U v {x) U T (x + 0) Ul(x + f) U$(x). Ther- 
malized lattice configurations can be generated by a com- 
bination of pure-SU(3) Cabibbo-Marinari updates fol- 
lowed by random G2 transformations, and overrelaxation 
sweeps, see p2j for a more detailed description. 

In most cases we work on asymmetric lattices with 
different lattice spacings in time (at) and space (a s ) di- 
rections, therefore the action (|3.7p contains a parameter 
Co called the bare anisotropy. To reproduce the usual 
classical gauge action in the naive continuum limit, the 
bare anisotropy has to be chosen as £ = a s /a t . However, 
what we need to achieve is a certain fixed ratio of lattice 
spacings in space and time directions in physical units, 
C = aP hys /a^ hys . The renormalized anisotropy C is a func- 
tion of Co and /3 in the Wilson action (|3.7p , below we will 
discuss a nonperturbative determination of Co = /(£>/?) 
to get a desired value of C for a given coupling (3. Before 
that, we will recall a few basic facts about the determi- 
nation of the static potential from rectangular Wilson 
loops. 



C. Static potentials from the lattice 

Let us start from a symmetric lattice with lattice spac- 
ing a. On a set of thermalized lattice configurations, we 
measure expectation values W(r, t) of Wilson loops of 
spatial extent r = r-a and a temporal separation t = t-a. 9 
For large time separations 

W(r, t) = c(r, a) exp[-V(r, a)t] for large t. (3.8) 



TABLE II: The bare anisotropy (,o(r,y) obtained from differ- 
ent (r, y) pairs entering Eq. (|3,17p . 



p 


Co (2, 2) 


Co (2, 3) 


Co (3, 2) 


Co (3, 3) 


9.5 


1.572(14) 


1.586(45) 


1.580(41) 


1.596(279) 


9.7 


1.626(22) 


1.638(47) 


1.627(50) 


1.631(155) 



The static potential 
V(r,a) 



lim - In W(r, t) 

t — >oo t 



(3.9) 



consists of two contributions: the true r-dependent static 
potential due to the interaction of a color source and 
anti-source, and the self energy contribution (indepen- 
dent of r) divergent in the continuum limit. Thus we can 
write 



V(r, a) = - V iat (r, a) + L sclt '(r, a) 
a 



(3.10) 



(Potentials with a hat are dimensionless.) 

On an asymmetric lattice one can introduce three dif- 
ferent potentials, defined as: 



V ss (r,a a ) 



V st (r,a s ) 



V ts (r, a t ) 



lim - \nW ss (x, y) 

y—fcc y 

— lim -\nW ss (x,y), (3.11) 

a s y^oo y 

— lim - In W s t (x, t) 

t — >oo t 

■— lim jhxW st (x,t), (3.12) 
at t—too t 

lim -\nW st {y,t) 

y—>oo y 

— lim - \nW st (y,t). (3.13) 

a s y~*oo y 



The Wilson loop W ss has both sides in spatial directions, 
while the second side of W st lies in the temporal direc- 
tion (i.e. the one in which the lattice is discretized more 
densely) . 

One can argue for a simple relation between the inter- 
action potentials [n| H3] : 



V;?(r,a s ) = V£\Zr,a t ). 
D. Renormalization of the anisotropy 



(3.14) 
(3.15) 



7 Another parametrization of G2 matrices, via Euler angl es |29| 
using real matrices, was described in Appendix B of Ref. I2ll . 

8 We are grateful to Axel Maas for deriving explicit expressions. 

9 Boldface symbols, x,t, . . ., denote variables with the dimension 
of length/time, x,t, . . ., are used for dimensionless quantities, i.e. 
x, t, . . ., divided by a, or a s , at on an asymmetric lattice. 



Arguments for the usefulness of working on an asym- 
metric lattice have been clearly discussed e.y. in floL f3ll| : 
a reduction of lattice artifacts, a certain hope to decrease 
errors due to the influence of the self energy V sclt : etc. 
The determination of static potentials at a fixed separa- 
tion of sources relies on relations like (|3.8p and requires 
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to evaluate W at as many t values as possible. Our main 
motivation for using asymmetric lattices thus was to have 
more data points available in the physical t window. We 
worked with lattices L 3 x (2L), but wanted to keep the to- 
tal physical length of the lattice in all directions the same, 



i.e. the renormalizcd anisotropy £ 



7 phys / a P h y s 



2. To 



determine the bare anisotropies £o = /(£,/?) that lead to 
£ = 2 at the used values of the coupling j3 we applied the 
procedure of [3l|. One requires that the static potentials, 
which one can introduce on an asymmetric lattice, ful- 
fill conditions (|3.14|) and (|3 . 1 5[) . To avoid the nontrivial 
task of reliable determination of AV in Eq. (|3 . 14[) . we 
estimated £o using the latter condition. 

In practice, we computed, for a given /3 and lattice size, 
quantities: 

p / \ W ss {r,y) W st (r,t) 
Rss (r, y) = , — r , R s t (r, t) =- 
W ss (r + l,y) 

and searched for a value of £o to fulfill 



w st (r + i,ty 

(3.16) 



R ss {r,y) = R st (r,t = 2y). 



(3.17) 



Of course, the obtained value of £o somewhat depends on 
the chosen pair of side lengths (r,y) in Eq. (|3.17[) . Since 
larger Wilson loops have smaller values and larger errors, 
we restricted ourselves to combinations (2 — 3, 2-^3) only. 

The described nonperturbative determination of £o was 
performed in simulations for j3 = 9.5 (about 3000 config- 
urations at each trial £o) and 9.7 (1000 configurations at 
each £o) on a 12 3 x 24 lattice. The results are summarized 
in Tabic |TTJ For production runs we chose values close to 
£o(2, 3), but checked a rather weak dependence of results 
on its precise value. The value chosen for (3 = 9.6 was 
obtained by linear interpolation. 

Our final set of bare anisotropies used in production 
runs was £ = 1.590,1.615,1.640, for (3 = 9.5,9.6,9.7, 
respectively. 



E. Smearing 

Another problem of the determination of static poten- 
tials from Wilson loops is the insufficient overlap of the 
trial quark-antiquark state with the ground state, espe- 
cially at larger spatial distances. An accepted remedy for 
this disease is not to construct Wilson loops directly from 
link matrices, but to use smeared links instead. Various 
variants of smearing have been proposed in the literature, 
we used the stout smearing procedure of Morningstar and 
Peardon [32|. The smearing is performed in several steps: 
First we calculate a weighted sum (x) of staples neigh- 
boring the link U^{x): 



S^x) 



P 



^[^(x)LV(x + i>)t/t(.T + A) 



p is a free parameter of smearing. We sum over di- 
rections v depending on the potential studied (perpen- 
dicular to the Wilson loop, used for calculation of the 
potential). In the second step we project the product 
= S,j.(x)U^(x) (no sum over p\) into the G2 alge- 
bra using: 



Qn(*) = J2 Im ( Tr [°^) *«])*«, 



(3.19) 



where t a are generators of the algebra (see Appendix [A"|l . 
In the last step, the smeared link is calculated as 



U s ™(x)=exp[iQ ll (x)}U^x) 



(3.20) 



The exponential function of matrices was calculated by 
means of spectral representation and LAPACK subrou- 
tines. 

After an extensive exploration of the effect of the free 
parameter p and of the number of smearing steps N sm 
used on the ground-state overlap, we chose p = 0.1. The 
optimal number of smearing steps was determined for 
each kind of potential: N sm = 20 for V s t, N sm = 15 for 
V ss , and N sm = 5 for Vt s (sec Appendix [Bl for details). 
In production runs we always applied all three smearing 
procedures on each thermalized configuration separately, 
i.e. each time we smeared a new configuration using 
different staples and different number of smearing steps 
to get best overlap for a given potential. 

As expected, smearing did enhance the ground-state 
overlaps considerably for all studied types of potentials. 
Not to overcrowd the main body of the paper with tech- 
nicalities, some illustrative examples are supplied in Ap- 
pendix [5] 



IV. G 2 ON A LATTICE: POTENTIALS 

Determination of potentials and setting the 
physical scale 



+ tf_ v {x)Up{x - P)U v {x -*> + £)!, (3.18) 



As discussed earlier fSec. HII Cj) . on an asymmetric lat- 
tice one can determine three kinds of potentials from the 
expectation values of Wilson loops of various orienta- 
tions, W\(r, t), where the index A assumes values ss, st, 
or ts. The corresponding potentials (in lattice units) can 
be determined [cf. Eq. (|3.8[) ] by fitting logarithms of the 
measured loops by linear functions in t: 

-]nW x {r,t) = C x + Vx(r)-t, A e {ss, st, ts}, (4.1) 

in an appropriate interval of t values, from t m j n to t max . 
Usually, to get a reasonable fit, one has to seek a compro- 
mise between two conflicting requirements: on one hand, 
to have a sufficiently long i-interval, i.e. to use as small 
i m in and as large £ max as possible, and, on the other hand, 
not to use too large t values to avoid loops measured with 
unacccptably large error bars. In practice, however, with 
smeared links the problem turned out to be less severe: 
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TABLE III: Simulation parameters. (Smearing parameters 
were listed in Sec. MI El In runs on symmetric lattices ther- 
malized configurations were separated by 20 sweeps consisting 
of 5 updates described in Sec. IIII Bl and 1 overrelaxation step, 
on asymmetric lattices we used configurations separated by 
15 such sweeps.) 
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FIG. 1: Logarithm of the expectation value of the fundamen- 
tal Wilson loop W s t(r, t) for fixed r- values as a function of t. 
(14 3 x 28 lattice, f3 = 9.6. Loops were computed from smeared 
links. Lines are drawn to guide the eye.) 



logarithms of Wilson loops were to a good approximation 
linear down to t = 1 (see below), so we could choose for 
our fits i m i n = 1. The upper limit t ma x was chosen as a 
value, at which the quantity used to monitor the ground- 
state overlap (see Appendix [B]) dropped considerably or 
its error was of the size of its value. 

A linear fit (|4.ip in a fixed (i m i n , i max )-interval using 
the standard least-square method provides an expected 
value of V\(r) together with an estimate of its statistical 
error, SV\(r). This is almost certain to be an under- 
estimate of the true error of the determination of the 
potential: there exists also a systematic error due to the 
unavoidably subjective nature of the choice of the fitting 
interval. To make at least a rough estimate of the order 
of magnitude of this systematic error, we used a kind of 
"jackknife-like" procedure: From the interval (1, i max ) we 
omit in turn the 1-st, 2-nd, etc., j-th point and make a 
linear fit (|4.ip to data points without the omitted point. 
In this way we obtain an ensemble of fit values for V\ (r) , 



TABLE IV: Parameters of the static potential between 
fundamental- representation charges, Eq. (|4.3p . symmetric lat- 
tice (14 4 ). The physical scale a B /ro was set using Eq. (|4.2|l . 
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FIG. 2: Potentials V,t, V ss , and Vts- Distances r are measured 
in units of a s . (14 3 x 28-lattice, /3 = 9.6, £ = 1.97.) 



and estimate the systematic error of the procedure as the 
square root of the variance of this ensemble. It turns out 
that at small r values the systematic error is of the order 
of the statistical one, but at large r values it is often an 
order of magnitude larger. 

Finally, if we want to compare potentials at different 
values of the coupling j3, we need to fix the physical scale. 
Usually, lattice spacings a s are expressed in units of the 
Sommer scale Tq defined via the relation [33| : 



,dV 
dr 



= 1.65. 



(4.2) 



r=r /a s 



In lattice QCD, one uses for the Sommer scale a phe- 
nomenological value of r = 0.5 fm, in the case of G2 we 
leave this scale unfixed. In fact, for illustrative purposes 
it suffices to use approximate proportionality of a s (/3) / rg 
to y/a((3), where t?(/3) is the string tension measured in 
lattice units (see Figs. [7] and [5] in Sec lIVCj) . A small cor- 
rection due to the Coulomb coupling constant does not 
change the picture qualitatively. 



B. Potential for fundamental representation 

In this subsection we summarize selected results for 
potentials between static quarks and antiquarks from the 
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TABLE V: Parameters of the static potential between fundamental-representation charges, Eq. (|4,3|l . asymmetric lattice (14 4 x 
28). The physical scale a 3 /ro was set using Eq. (|4.2[) from V s t- 
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FIG. 3: Differences between potentials V sa and Vst (multiplied 
by £). (14 3 x 28-lattice: (3 = 9.5, £ = 2.00; = 9.7, £ = 1-96.) 

fundamental representation, and the main purpose is to 
calibrate our procedure and compare with results of an 
earlier investigation [2l|. Results for higher representa- 
tions are deferred until Sec. IIVC1 

Our simulations were performed at three values of the 
coupling j3 on the weak coupling side of the crossover re- 
gion observed in Ref . [l7| . Parameters of the production 
runs are summarized in Table HTT1 

Fig. [1] clearly illustrates the fact that has been men- 
tioned in Sec. II V Al logarithms of Wilson loops computed 
from smeared links are almost perfect linear functions 
of t. The example shows W s t(r,t) for a series of t val- 
ues at fixed r. The corresponding value of the potential 
was then obtained by a fit (|4.1[) . The resulting potentials 
were then parametrized by the usual Coulomb plus linear 
form: 

V x (r) = c x +d x r, \£{ss,st,ts}. (4.3) 

r 

(A hat over a number again indicates that it is dimcn- 
sionlcss lattice quantity, no hat is needed on a.) 

To compare with results in the existing literature, we 
computed a and a in runs on a symmetric, 14 4 lattice. 
The result is summarized in Table IIVI The values of the 
strin g te nsion are in reasonable agreement with those of 
Ref. [21(, obtained on 12 4 lattice. The slight discrepancy 
between a values can be attributed to different smear- 
ing procedures adopted. (Ref. [2l[ used a variant of the 



procedure of [3J|.) The purpose of smearing is to reduce 
ultraviolet noise in lattice configurations that mostly in- 
fluences the self-energy and Coulombic parts of the static 
potential, and efficiency of various procedures differs in 
achieving this goal. At (3 = 9.7 also string tensions seem 
not to be compatible within errors. However, both num- 
bers are quoted with statistical errors only, and may in 
fact agree if finite-size effects and/or systematic errors 
are taken into account. To support this claim, we deter- 
mined the string tension by our procedure also on 12 4 
lattice, and obtained a value of 0.109(2). Finally, using 
irnin = 2 instead of 1 in potential fits (c/. the discussion 
in Sec. IIV A[) . we went further down to a = 0.105(2). It 
should be clear that though statistical errors on string 
tensions are as low as 2-3%, the actual error may be of 
the order of 10%. 

Next we turn to asymmetric lattices and determine 
parameters of the three potentials (|3.1ip . (|3 . 1 2|) . and 
(|3.13j) . Those are summarized in Table Ivl 10 When look- 
ing at these numbers, one should keep in mind that the 
errors quoted are just statistical. We did not attempt 
to estimate systematic errors, which were discussed in 
Sec. IIV Al The most reliable determination of parame- 
ters of the static potential comes from Wilson loops W s t 
with the i-direction of the loop along the long lattice di- 
rection (in which the lattice is more finely discretized). 

One can now verify the correctness of the chosen values 
for the bare anisotropy £o (c/. Table HI)) . On the basis of 
Eq. (|3.14p one expects a ss ~ 2a s t and a ss ~ 2a s t, which 
is satisfied to a surprising accuracy. From the numeri- 
cal values of the string tensions one obtains the renor- 
malized anisotropies £ = 2.00(7), 1.97(12), and 1.96(16) 
for j3 = 9.5, 9.6, and 9.7, respectively. Similar relations 
should hold between ss and ts parameters [sec Eq. (13.151) ]. 
However, they are fulfilled with worse accuracy, because 
of systematic errors mentioned above. The available t 
intervals for determination of these potentials are short, 
since the overlap with the ground state drops fast with t 



We have some data on static potentials for the fundamental rep- 
resentation also from simulations on a 16 3 X 32 lattice, albeit 
with lower statistics. The estimated string tensions agree with 
those quoted above, no strong finite-size effect was observed in 
this case. E.g., our preliminary values from the 16 3 X 32 lattice 
are a st = 0.199(2), 0.129(2), 0.097(2) at /3 = 9.5,9.6,9.7, respec- 
tively. 
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at larger values of r (see Fig. [TT] in Appendix [B| . 

The fulfillment of relations (|3.14| and (|3.15|) is visually 
displayed in Fig. [5] for (3 = 9.6. The difference between 
V ss and £ V s t is expected to be a constant independent 
of r (difference between self energies) , our data support 
the expectation, see Fig. [3] 

Finally, we checked that the exact value of the bare 
anisotropy £o is not of vital importance on asymmetric 
lattices. There were no large differences between the ob- 
tained potentials and their string tensions when we used 
£o = 1-61 instead of 1.59 at (3 = 9.5, or £o = 1-62 instead 
of 1.64 at = 9.7. 



C. Higher representations 

After subjecting our procedure for determination of 
potentials to a number of tests in the case of the fun- 
damental representation, we are ready to proceed to the 
main subject of the present paper: potentials between 
static color sources from higher representations of the 
G2 group. For this purpose we will concentrate on po- 
tentials V s t which can be determined most accurately. 
The procedure is now in principle straightforward. A 
thermalized configuration on an L 3 x (2L) lattice is a 
set of link matrices in the fundamental representation. 
We smeared these links over "short" directions using the 
procedure described in Sec. IIII El Then, we compute a 
fundamental Wilson loop as a product of smeared link 
matrices along a closed loop (one side in a "short" direc- 
tion, the other in the "long" one). Before taking a trace, 
we compute the matrix representing this loop in the ad- 
joint representation using Eq. (|3.4[) . Then we trace both 
the fundamental and adjoint Wilson loop matrix, and 
from these traces determine also those of the loop ma- 
trices in higher representations via the set of formulas 
()3.3p . Finally, to obtain the expectation values of traces 
of Wilson loops of size (r x t), we average over loop po- 
sitions on a lattice, and over an ensemble of thermalized 
configurations in the usual way. 

Illustrative results for Wilson loops in higher represen- 
tations are displayed in Fig. 2] for (3 = 9.6. The results 
look very much the same as in the case of the funda- 
mental representation (Fig. [T|). Values of the potential 
at various separations r are obtained by fits of the form 
(14. ip . The results for (3 = 9.6 are shown in Fig. [5l It is 
clear that all higher potentials show a linear rise similar 
to the fundamental representation in the same range of 
distances. 11 



11 Though the string breaks in all representations of G2, we had no 
ambition to observe it in our simulations. It is well known from 
earlier studies (see e.g. |35|Q that it is of utmost difficulty to see 
string breaking on a lattice. Such a task is currently far beyond 
our possibilities, and is also not relevant to the main message of 
the paper. 



The obtained potentials were fitted again by the 
Coulomb plus linear form (|4.3|) . The resulting values of 
string tensions are summarized in Table IVII The mea- 
sured values lie very close to predictions based on Casimir 
scaling, and the agreement with Casimir scaling tends to 
improve by increasing (3, i.e. the deeper we enter the 
weak-coupling region above the crossover. In Fig. [5] this 
fact is exemplified on ratios of different representation 
potentials to the fundamental one at fixed r. At all dis- 
tances, the ratios are close to the Casimir-scaling predic- 
tion, even though we made no attempt to get rid of the 
self-energy contributions. 

In this context, the question of systematic errors again 
pops up: Could they not invalidate the above agreement 
with Casimir scaling? As discussed in detail in Sec. lIV Al 
the main source of the systematic error is the choice of 
the t- interval for linear fits of the form (|4.1[) . To produce 
results in Table PVTI t m - m = 1 was used. The same analysis 
with t m i n = 2 for the three lowest representations leads 
to values given in Table IVII1 The Casimir-scaling result 
comes out to be quite robust against this modification: 
though string tensions changed in all representations by 
up to 10%, they went in the same direction, and thus 
their ratios changed only slightly. In fact, at the largest 
available /3, the agreement with Casimir-scaling predic- 
tion even improved. 

Finally, we display the static potential for color charges 
in the fundamental (Fig. [7]) and adjoint representation 
(Fig. [5]) together for all three (3 values, expressed in phys- 
ical units defined by the fundamental string tension. All 
potentials nicely follow a universal curve. 



V. CONCLUSIONS 

Alice, in Lewis Carroll's Through the Looking- Glass, 
and What Alice Found There, enters a garden, where 
flowers speak to her and mistake her for a flower. Does 
the gauge theory with the exceptional group G2 belong 
to the same species of flowers with "ordinary" confin- 
ing SU(iV) gauge theories, or is it different and only 
mistaken for a flower? Many results of recent investi- 
gations suggest the answer that all non-Abelian gauge 
theories, including G2, share many properties in the in- 
frared. Thermodynamic properties are similar and 
Ref. [22| provided further numerical evidence: a quali- 
tative agreement among Landau-gauge ghost and gluon 
propagators and the Faddeev-Popov operator eigenspec- 
trum in SU(2), SU(3) and G2 gauge theories, though only 
in two and three Euclidean dimensions. 

From the point of view of the center-vortex confine- 
ment mechanism, theories also behave similarly. When 
center vortices are present and percolate, in theories with 
non-trivial center, they cause the expected iV-ality de- 
pendence of asymptotic string tensions for static poten- 
tials for color charges from different representations of the 
gauge group. In a theory with trivial center, the asymp- 
totic string tension for all representations is zero. But 
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FIG. 4: Logarithm of the expectation value of the Wilson loop Wst(r,t) for fixed r- values as a function of t, in the adjoint, 
{27}-, {64}-, and {77}-representations. (14 3 x 28 lattice, f3 = 9.6. Lines are drawn to guide the eye.) 



TABLE VI: String tensions a a t for different representations. (14 3 x 28 lattice. Because of large errorbars, it was not possible 
to determine reliably string tensions for the {77} and {77'} representations at (3 — 9.5.) 
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even if a centerless theory possesses no non-trivial center- 
vortex-like objects, its vacuum can contain domain struc- 
tures. Ref. [21( suggested a universal and simple model in 
which such vacuum domains - assuming color magnetic 
fields fluctuating almost independently in each domain, 
only fulfilling a constraint that the total flux corresponds 



to an element of the center - cause Casimir scaling of 
static potentials in the intermediate distance interval, un- 
til color screening sets on. Casimir scaling was observed 
in lattice simulations of SU(2) and SU(3) gauge theories 
long ago, the model prediction called for verification also 
in the G2 lattice gauge theory. 
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TABLE VII: Same as Table EH from potential fits (gTj with t min = 2. 
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FIG. 5: The potentials V s t for different representations vs. 
the dimensionless r. (14 3 x 28 lattice, /3 = 9.6.) 
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APPENDIX A: G 2 GENERATORS 

The G2 algebra for our representation (based on p8j ) 
is defined by the following generators: 



Our results convincingly demonstrate (approximate) 
Casimir scaling for static potentials between color 
charges from the six lowest representations of G2. 
The agreement between measured values of intermedi- 
ate string tensions with predictions based on values of 
quadratic Casimirs is quite striking, they differ by at 
most 10-15%, and this can hardly be just a numerical 
coincidence. The results of course cannot prove that the 
model is right, but combined with the solid evidence for 
Casimir scaling in SU(2) and SU(3) provide support for 
its main ingredients (common for all groups) - magnet- 
ically disordered vacuum with a domain structure. How 
to identify such domains is another question; in theo- 
ries with non-trivial center one can apply the method of 
center projection in the maximal center gauge (36[, thus 
fitting extended ( "thick" ) vacuum structures by thin cen- 
ter vortices. The identification of vacuum-type domains 
with trivial total flux remains problematic. 



The observed Casimir scaling, coming out so naturally 
within models with magnetic disorder and vacuum do- 
main structure, may pose a challenge to approaches that 
invoke disparate ideas to explain the phenomenon of con- 
finement. 
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FIG. 7: The fundamental-representation potential in physical 
units, all /? values, from 14 3 x 28 lattice. (The constant c, cf. 
Eq. (|4.3[) . is subtracted. The line is drawn to guide the eye.) 
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FIG. 6: Ratios V st ,{D} /V s t,F for different representations {D} 
are shown as functions of the dimensionless r. Horizontal 
lines show Casimir-scaling predictions. (14 3 x 28 lattice. A 
few points from off-axis Wilson loops are also included, those 
were not used in string-tension determination.) 
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APPENDIX B: OVERLAP WITH THE GROUND 
STATE 

In optimizing our smearing procedure and monitor- 
ing the overlap with the ground state, we were inspired 
by [I?]] • We investigate, at fixed r, a function of expec- 
tation values of Wilson loops defined as: 
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(Bl) 



This is an approximant to the ground-state overlap (to be 
reached in the t — > 00 limit), and decreases monotonically 
(in t) to its asymptotic value. So an indication of the 
closeness to asymptotia is when the quantity reaches a 
plateau as a function of t. 

Another quantity of interest is: 
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which approaches, from above, the ground-state potential 
for t — > 00. 

For an asymmetric lattice, we can define three different 
types of quantities Co(r, t), corresponding to three differ- 
ent potentials: 
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lattice, /3 = 9.6.) 
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FIG. 11: The quantity co(r, t) as a function of r, in the smeared 
case for the fundamental representation. (14 3 x 28 lattice, 
= 9.6.) 
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FIG. 12: The quantity co,at{r,t) as a function of t in the 
smeared case for the representations {14}, {27}, {64}. (14 3 x 
28 lattice, j3 = 9.6.) 



r= 1 
r = 2 
r = 3 



and likewise three types of M e s(r, t). 

The smearing procedure described in Sec. IHI El de- 
pends on the parameter p and the number of smearing 
steps N sm . Our aim was to tune the set to maximize the 
ground-state overlap and to minimize M e g. Typically, 
the optimal number of smearing steps grows by decreas- 
ing the value of p. After a few trials we fixed p = 0.1 and 
looked for the optimal value of N sm . The typical behav- 
ior of CQ, st (r,t) and M e g tBt (r,t) for three sets of (r, t) is 
displayed in Fig. [5J In this particular case, the optimal 
number of smearing steps comes out to be N sm ~ 20, and 



depends only weakly on /3 and £o- In a similar way we 
found the optimal values of iV sm quoted in Sec. lIIIEl 12 



12 One should keep in mind - quoting from an extensive study of 
smearing by Diirr [3H — that "any smearing recipe drastically 
reduces the noise and a detailed tuning of parameters and/or 
iteration level is neither needed nor useful. [...] Only at absurdly 
high levels systematic effects get visible, presumably due to the 
reduced locality in terms of the original links." 
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As a final illustration, we present here the behavior 
of the achieved ground-state overlaps as functions of t 
for (3 = 9.6. (For (3 = 9.5 and (3 — 9.7 we observed 
the same behavior, only errorbars were slightly differ- 
ent.) Results are shown in Figs. [101 E] and [T2j In all 
figures only data with reasonably small errors are shown, 
for the sake of readability. In Fig. [TU] and QT] the ap- 
proximants of ground-state overlaps are shown in the 



unsmeared and smeared cases, respectively, for the fun- 
damental representation: we see clearly that the smear- 
ing procedure significantly increases the overlaps and the 
length of plateaus for all potentials. In Fig. [T^lwe display 
the behavior of co, s t{r, t) for higher representations and 
see a similarly good enhancement of the overlap with the 
ground state. Of course, the overlaps worsen when going 
to higher representations. 
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